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We investigate how imposing kinetic restrictions on quantum particles that would otherwise hop freely on a 
two-dimensional lattice can lead to topologically ordered states. The kinetically constrained models introduced 
here are derived as a generalization of strongly interacting particle systems in which hoppings are given by 
flux-lattice Hamiltonians and may be relevant to optically driven cold-atom systems. After introducing a broad 
class of models, we focus on particular realizations and show numerically that they exhibit topological order, 
as witnessed by topological ground-state degeneracies and the quantization of corresponding invariants. These 
results demonstrate that the correlations responsible for fractional quantum Hall states in lattices can arise in 
models involving terms other than density-density interactions. 


I. INTRODUCTION 

The pursuit of topological states of matter has fueled sub¬ 
stantial experimental and theoretical developments in physics. 
From the integer quantum Hall effect 1 to topological insu¬ 
lators 2,3 and beyond, topological states possess remarkable 
properties of both fundamental and technological interest. For 
instance, certain correlated topologically ordered states, such 
as fractional quantum Hall (FQH) states 4 , have the poten¬ 
tial to be used for fault-tolerant quantum computation''. The 
functionality required for novel high-end technological appli¬ 
cations, however, steers material design towards microscopic 
structures with increasingly complex features. 

Several steps have been taken in order to bridge the gap 
between theoretical descriptions of FQH-type topological or¬ 
der and experimentally accessible systems. Following the first 
theoretical proposal suggesting the presence of FQH states in 
a lattice model 6 , FQH physics was predicted to arise in the 
context of optically trapped cold atoms 7-9 . More recently, a 
class of lattice models called Chern insulators 10 have been 
shown to develop FQH-like ground states, termed fractional 
Chern insulators (FCI) 11-13 , upon introduction of a short- 
range repulsion between particles occupying a fractionally 
filled band, even in the absence of an external magnetic field 
(see Refs. 14,15 for comprehensive reviews). Initially, these 
models were fine tuned so that the lowest band of their energy 
dispersion imitates the lowest Landau level: It is almost per¬ 
fectly flat and topological, i.e., characterized by a Chern num¬ 
ber C = ±1 and hence named Chern band. Proposals for re¬ 
alizations of FCI states based on oxide heterostructures 16 , lay¬ 
ered multiorbital systems 17,18 , optical lattices 19 , and strained 
or irradiated graphene 20,21 followed soon thereafter. 

As bands in solids are generally not flat and interaction 
strengths vary, an obvious way to bring FQH-like states closer 
to reality is to relax the energetic conditions imposed on rele¬ 
vant systems in order to emulate Landau levels. It was rec¬ 
ognized early on that interaction strengths may be allowed 
to increase beyond band gaps 12,17,22,23 , and subsequent results 
showed that a finite dispersion may actually favor certain FCI 
states 18,24 . One can then venture into the strong-correlation 
regime by allowing for arbitrarily strong repulsion strengths. 
Surprisingly, one still finds robust FCI states — as well as 
more exotic, topologically ordered states 25 — even though in¬ 


teractions now mix bands with opposite Chern numbers 12,26 . 
Strongly correlated materials without sharply defined bands 
may therefore be considered as candidates for the realization 
of FCI or similar states. 

When particles repel their neighbors very strongly, it is rea¬ 
sonable to approximate the interaction as a hardcore constraint 
which does not allow particles to occupy neighboring sites. 
On two-dimensional lattices of corner- or side-sharing trian¬ 
gles, the configurations allowed by this constraint are identical 
to those of the so-called hard-hexagon (HH) model of classical 
lattice gases 27,28 . In the quantum version of these HH models, 
transitions from one configuration to another are caused by 
hopping terms. When the hoppings considered give rise to 
Chern bands, imposing a hardcore constraint in partially filled 
lattices can lead to FQH-like ground states 26 . 

Here we explore the possibility of obtaining the same 
physics by replacing the hardcore constraint with a kinetic 
one. Instead of energetically penalizing or — in the infinite- 
interaction limit — removing configurations from the Fock 
space altogether, we start with all configurations being a priori 
equivalent. Constraints are then introduced only in the tran¬ 
sitions between configurations. The resulting systems can be 
thought of as quantum versions of cooperative lattice gases 29 . 
Similarly constrained quantum particles have been studied in 
several contexts in the past, with the kinetic constraints often 
giving rise to nontrivial correlation effects 30-33 . In this paper, 
we show that when the constrained kinetic terms are endowed 
with appropriate Berry phases that break time-reversal sym¬ 
metry, the obtained ground states show definitive features of 
FQH topological order. 

This result is far from being purely academic: Cur¬ 
rent experiments on optical lattices can generate a partially 
filled Chern band in the laboratory, with tunable interaction 
strengths 34,35 . Some of the terms that arise in the model¬ 
ing of the experimental setups exhibit precisely the type of 
constraints studied here 34,36 . The same kinetic constraints are 
also an approximate limit of strong short-range repulsion and 
therefore relevant to the realization of FQH physics in strongly 
correlated materials, as will be shown below. 

The rest of this paper is structured as follows: In Sec. II we 
present a general formalism of tightly-bound particles inter¬ 
acting via strong nearest-neighbor repulsion. In the infinite- 
interaction limit particles obey a hardcore nearest-neighbor 
exclusion principle, motivating us to re-interpret the prob- 
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lem as a quantum hard-hexagon model. We then relax the 
static exclusion constraint to a kinetic one, derive two exam¬ 
ples of constrained fermionic models, and pose the question: 
Can topological order arise from the kinetic constraints? In 
Sec. Ill we investigate the topological nature of the ground 
states of the constrained models using exact diagonalization. 
The results we obtain prompt us to answer the question affir¬ 
matively. Our conclusions are summarized in Sec. IV. 


II. MODELS: FROM STRONG INTERACTIONS TO 
KINETIC CONSTRAINTS 

A. General formalism 

To establish a general formalism, we consider N quantum 
particles of a single species that hop on a two-dimensional 
lattice A. The discussion below will be limited to spin¬ 
less fermions, but there is no fundamental obstacle to treat¬ 
ing bosons in the same fashion. Distances between nearest- 
neighboring sites of A are set to unity for convenience, and 
periodic boundary conditions are assumed in both spatial di¬ 
rections. The general form of the Hamiltonians we consider 
below reads 

^ = E(^4s+ h - c -) . (d 

where c^ct) is a regular particle annihilation (creation) op¬ 
erator acting at position i of A. An appropriate choice of 
the operator-valued function AT can result in any possible 

n-body intersite term. For example, setting b - = 1 yields 
all possible hopping terms with equal amplitude, whereas 
4 = c i4 y> elds ad possible density-density interactions 
with equal magnitude. 

In the following, we restrict the F .T to products of hole- 
density operators of the form (1 — n) acting in the neighbor¬ 
hood of i and j. Such terms emerge effectively from strong 
interactions 37 . To see this, consider particles that hop around 
the lattice and at the same time interact with each other via 
a nearest-neighbor repulsion of strength V. In the infinite¬ 
ly limit, particles cannot occupy nearest-neighbor sites. This 
suggests that states containing nearest-neighboring particles 
can be removed from the Fock space. Accordingly, for any 
site A we define the projected operator c\ by demanding 
that its action on any state is to create a particle at i if and only 
if this site and all of its nearest neighbors are empty. Formally, 

4 ■■= 4 n (1 - hj) , (2) 

b' : lj-*|=i} 

where hj := ct Cj. If we allow particles to hop while obey¬ 
ing the above hard-core condition, then the system can be de¬ 
scribed by Eq. (1) with 

4 c = ^ n a - Ai) > (3) 

{Z:|J-i|=lV|I-j|=l} 


where t- are generally complex-valued hopping amplitudes. 
The strong-repulsion limit can therefore be recast into a 
Hamiltonian of the form of Eq. (1), with AT being a product 
of hole-density operators in the vicinity of i and j. On two- 
dimensional lattices of corner- or edge-sharing triangles, the 
allowed states are exactly the configurations of classical HH 
models, widely studied in the context of glassiness. Due to 
this resemblance, we shall call the above Hamiltonians quan¬ 
tum hard-hexagon (QHH) models. 

We now wish to reduce the above hardcore constraint to a 
new one that does not correspond to a density-density inter¬ 
action. One motivation for doing so is that density-hopping 
terms of the form m(l — n)c have been shown to arise in 
the description of optical lattice experiments, where trapped 
atoms are periodically driven using circularly polarized light 
or equivalent settings 34,35 . The theoretical treatment of rele¬ 
vant models describes how the driving, apart from affecting 
the preexisting hopping and density-density terms, introduces 
new frequency-dependent density-hopping terms 36 . The rel¬ 
evance of these terms in the experimental setups depends on 
the interaction strength between particles, which is in princi¬ 
ple tunable. Nevertheless, their precise effect in that setting 
requires careful study, as it may very well be scrambled by 
other processes. As a first step, however, it is of intrinsic inter¬ 
est to isolate the kinetically constrained terms and determine 
their properties in a simpler context. As we shall show below, 
such terms can generate nontrivial behavior even on their own. 

Notice that f ' HC is just a product of hole-density operators 
at different sites. In any partially filled system, products of 
(1 — n) operators that act on different sites are more likely 
to vanish than single (1 — n) operators. We can therefore at¬ 
tempt to restrict the product in Eq. (3) to only a few (1 — fi) 
operators and see whether this captures the same physics as 
the full product. To decide which terms to truncate, we draw 
inspiration from classical models of diffusion. As mentioned 
above, F HC is the quantum counterpart of the HH model on 
the triangular lattice. In the HH model, each particle can be 
visualized as a hard disk of radius 1/2 < r < 1, so that con¬ 
figurations with particles on nearest-neighboring sites are for¬ 
bidden. If the disk radius is reduced to v/3/4 < r < 1/2, then 
all particle configurations are a priori allowed. The particle 
motion, however, is still constrained if we assume that hard 
discs move from site to site along straight lines: On lattices 
of edge- or corner-sharing equilateral triangles, they cannot 
move past one another (see Fig. 1). This is an example of a 
cooperative lattice-gas model 38 . 

A quantum analog of the kinetically constrained model de¬ 
scribed above can be straightforwardly constructed by dis¬ 
carding all factors in Eq. (3) apart from those which pertain 
to common neighbors of i and j. With this choice, 

4 c =% n (! - «*) - (4) 

{Z:|i—i| = |«-j|=l} 

Examples of the resulting terms are pictorially represented in 
Fig. 1 for two lattices. Following the classical terminology, we 
shall call the processes represented by F^ c vacancy-assisted 
hoppings (VAH). Note that this type of constraint is not spe¬ 
cific to lattices containing triangular plaquettes: Next-nearest 
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FIG. 1. (Color online) Illustration of vacancy-assisted hopping mod¬ 
els on (a) the triangular lattice, with nearest and third-nearest neigh¬ 
bor hoppings as defined in Refs. 17 and 18, and (b) the kagome lat¬ 
tice, with nearest and second-nearest neighbor hoppings as defined 
in Ref. 39. Solid and dashed circles denote particles and vacancies, 
respectively. The hoppings (arrows) are only allowed if a vacancy is 
located at the position shown. Note that some of the hoppings are 
imaginary; see Eqs. (5) and (6). 


neighbor hoppings on checkerboard or honeycomb lattices 
can be similarly constrained. Here we focus on the triangular 
and kagome lattices because all hoppings present in the mod¬ 
els introduced below are constrained in the manner enforced 
by Eq. (4), and also because these models are known to yield 
the most robust fermionic FCI states. 


B. Example models 
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FIG. 2. (Color online) Flow of eigenvalues under insertion of mag¬ 
netic flux corresponding to a phase ip 2 through one of the handles of 
the toroidal system for (a),(b) the triangular and (c),(d) the kagome 
lattice, (a), (c) correspond to QHH models (i.e., F = F HC ) and (b), 
(d) to VAH models (i.e., F = F KC ). 


We now ask whether the correlations induced by the 
vacancy-assisted hoppings described above can generate topo¬ 
logical order. To answer this, we shall impose this kinetic 
constraint to tight-binding models with topologically nontriv¬ 
ial bands. We choose t- that correspond to two thoroughly 
studied FCI models, diagonalize the constrained Hamiltonian 
exactly on finite periodic clusters and look for signatures of 
topologically nontrivial states. 

Tight-binding Hamiltonians are represented as H := 

E k V’ltffcV'fc, where A = ( c i,fc. c 2 , fc . •••> c„, fc ) T is the spinor 
of annihilation operators on each of the n sublattices of the 
lattice model at momentum k. The triangular-lattice model 
introduced in Refs. 17 and 18 is given by 

H k - 9 k r, (5a) 

3 

90, k = 2i 3 cos ( 2fc ' a i ). ( 5b ) 

3 =1 

91, k = 2 *cos(fe • aj, * = 1,2,3, (5c) 

where a 1 = (1/2, —x/3/2) T , a 2 = (1/2, v/3/2) T , a 3 = 
— (a 1 + a 2 ) and t = (r 0 , r 3 , t 2 , t 3 ) is the vector of Pauli 
matrices including the 2x2 unit matrix as r 0 . The kagome- 
lattice model introduced in Ref. 39 is defined as 

H*~d k -\ + d k -\, (6a) 

d ik = —2t cos(fe • a,) + 2* 2 cos(fe ■ 6J, * = 1,2,3, 

(6b) 


d ik = — 2d cos (fe • a,) + 2(5 2 cos(fc • b/, * = 1,2,3, 

(6c) 

where a 1 = (1,0) T , a 2 = (1/2, -\/3/2) T , a 3 = a 2 — a 3 , 
b 1 = a 2 + a 3 , b 2 = a 3 — a 3 , b 3 = a 1 + a 2 and A = 
(AJJA 2 JA 3 ), A = (A 1 ,A 2 ,A 3 ) are vectors of the Gell-Mann 
matrices 

/0 1 0 \ /00 l \ /0 0 0 \ 

A 1 =(l 0 0 I , A 2 = [ 0 0 0 I , A 3 = [ 0 0 1 | , 

\o o o/ V 1 0 °/ \° 1 °/ 

(7a) 

/ 0 i 0\ /0 0 — i\ /0 0 0\ 

A 1= -i 0 0 ,A 2 = 0 0 0 ,A 3 = 0 0 i . 

\ 0 0 0/ \i 0 0 / \0-i0/ 

(7b) 

We set t 3 /t = 0.19 for the triangular lattice and t 2 /t = 
0.3, S/t = 0.28,(5 2 /* = 0.2 for the kagome lattice. These val¬ 
ues of the hopping terms are chosen to yield relatively flat low¬ 
est bands with Chern number C = — 1. We have verified that 
all our results remain valid in a finite range around the values 
chosen here. We focus on average particle densities p = 1/6 
for the triangular lattice and p = 1/9 for the kagome lattice, at 
which the lowest band of the pure hopping counterparts of the 
models we study here would be at filling v = 1/3. We then 
impose the constraint encapsulated in Eq. (4). The resulting 
vacancy-assisted hopping models are sketched in Fig. 1. 
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FIG. 3. (Color online) (a) Berry curvature of the triangular-lattice 
model on a 42-site cluster with N = 7 as a function of the mag¬ 
netic fluxes ipi and <p 2 threading the two handles of the toroidal sys¬ 
tem; (b) error in the quantization of a H as a function of the number 
of subintervals in the partition of the Brillouin zone of fluxes. The 
unidirectionality of the stripes along ip i is merely an artifact of the 
choice of the alignment of the two-site unit cell of the model [see 
Fig. 1(a)]. 


FIG. 4. (Color online) Gap to excited states at zero flux as a func¬ 
tion of inverse number of lattice sites N g for QHH (solid triangles) 
and VAH (empty triangles) on (a) the triangular-lattice model and (b) 
the kagome-lattice model. The value of the gap is the energy differ¬ 
ence between the highest in energy quasidegenerate FCI level and 
the lowest level that does not belong to the FCI state manifold. 


HI. RESULTS 

FQH-like topological order at v = 1/3 manifests itself 
in the many-body energy spectrum as spectral flow, which 
is the exchange of the threefold quasidegenerate many-body 
ground-state levels upon insertion of a flux quantum through 
one of the handles of the toroidal (i.e., periodic) system. This 
spectral property is shown in Fig. 2. To establish consistency 
with previous results 26 , Fig. 2(a) shows the spectral flow of 
the QHH version of the triangular-lattice model, where simul¬ 
taneous occupation of nearest-neighboring sites is prohibited 
but hoppings are unconstrained. In Fig. 2(c) we show that 
the same physics occurs in the kagome-lattice model. Apart 
from the characteristic FQH features in the energy spectra, the 
ground states obtained have a very precisely quantized Hall 
conductivity <j h = e 2 /(3 h). 

We now reduce the HH constraint of Eq. (3) to the kinetic 
constraint of Eq. (4), as outlined above. We notice that the 
energy spectra, shown on the right of Fig. 2, are qualitatively 
equivalent. The symmetry sectors in which the quasidegener¬ 
ate ground states reside are the same as in the case of hard¬ 
core interactions and can be predicted by methods already in 
use for finite FQH and FCI systems 13,40 . More importantly, 
the ground states of the VAH models have nontrivial topo¬ 
logical characteristics. Their many-body Berry curvature is a 
smooth function that integrates to a precisely quantized Hall 
conductivity, even for the finite systems considered here. An 
example of this quantization is presented in Fig. 3 for the 
triangular-lattice model. Equivalent results can be obtained 
for the kagome-lattice model as well. 

The FQH-like ground states of the VAH models are generi- 
cally less gapped than those of their QHH counterparts, as can 
be seen in Fig. 4. Nevertheless, there is a finite volume in pa¬ 
rameter space in which one obtains gapped FQH-like ground 
states even in the VAH models for most clusters. Finite-size 
effects — such as the noticeable energy separation between 
quasidegenerate ground-state levels — are sizable and do not 


allow for a definitive answer as to whether there is topological 
order in the ground state of the VAH models in the thermody¬ 
namic limit. We note, however, that the gaps of the VAH mod¬ 
els do seem to follow the same trend as those of their QHH 
counterparts, which remain well gapped up to the largest ac¬ 
cessible system sizes (see also Ref. 26). The finite-size effects 
observed in our calculations are similar to the ones observed in 
previous, more detailed studies of the dependence of FCI state 
energies on cluster size and aspect ratio, which concluded that 
these effects are indeed irrelevant to the underlying physics 41 . 
Furthermore, we find no signatures of a charge instability, so 
the only evident competitor for the ground state is a compress¬ 
ible metallic state. The latter is clearly disfavored according 
to previous detailed studies of FCIs in the models defined by 
Eqs. (5) and (6) in the parameter regimes chosen here 18,42 . 


IV. CONCLUSION 

Regardless of the detailed energetics of the models derived 
above, the key conclusion of this work is that kinetic con¬ 
straints in partially filled Chern bands can generate topolog¬ 
ically ordered states. This can be intuitively understood by 
noting that kinetic constraints can be seen as an approximate 
generalization of hardcore nearest-neighbor repulsion. The 
ground states of the latter can be adiabatically connected to 
FCI states, which in turn means that correlations crucial for 
FQH-like states can be generated by kinetic constraints alone. 
This may be a stepping stone for further intuition into the mi¬ 
croscopies of the FQH effect itself. More hints in this direc¬ 
tion can be found in recent work 43 , where the intimate relation 
between FQH and FCI states gives rise to similar density¬ 
hopping terms. Finally, we briefly comment on choosing to 
introduce fermionic models in this work. The reason for our 
choice is that the v = 1/3 fermionic FCI state is the most 
robust state — and therefore the easiest one to study — that 
requires short-range interactions. If one is interested in mod¬ 
eling current experiments precisely, it would be more mean- 
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ingful to pursue bosonic FCI states. However, the most robust 
bosonic FCI state at v = 1/2 requires only an on-site inter¬ 
action, whereas the v = 1/4 state requires longer-range re¬ 
pulsion, which would consequently give rise to less intuitive 
kinetic constraints. Since here we are interested in whether ki¬ 
netic constraints can in principle induce topological order, we 
have focused on the simpler case of the v = 1/3 fermionic 
FCI. Nevertheless, similar tendencies towards topological or¬ 
dering are to be expected when suitable kinetic constraints are 
introduced in bosonic versions of Haldane-like models 44 , with 
potential repercussions on the realization of topologically or¬ 


dered cold-atom systems 34 36 . 
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